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Abstract 

In this paper, we study a nonlocal degenerate parabolic equation of order a + 2 for 
a £ (0, 2). The equation is a generalization of the one arising in the modeling of hydraulic 
fractures studied by Imbert and Mcllet in 2011. Using the same approach, we prove the 
existence of solutions for this equation for 0 < a < 2 and for nonnegative initial data 
satisfying appropriate assumptions. The main difference is the compactness results due to 
different Sobolev embeddings. Furthermore, for a > 1, we construct a nonnegative solution 
for nonnegative initial data under weaker assumptions. 


1 Introduction 


this paper, we study the following problem 



1 

f d t u + d x (u n d x I(u)) = 0 

for x € f2, 

t > 0, 


d x u = 0,u n d x I(u ) = 0 

for x £ dfl, 

t > 0, 

1 

[u(0, x) = u 0 (x) 

for x £ f2, 



where C = (a, b) is a bounded interval in R, n is a positive real number and I is a nonlocal elliptic 
negative operator of order a defined as the a/2 power of the Laplace operator with Neumann 
boundary conditions I = —(—A)t where a € (0, 2); this operator will be defined below by using 
the spectral decomposition of the Laplacian. 

The case a = 1 was studied by Imbert and Mellet [15] who proved the existence of non¬ 
negative solutions for nonnegative initial data with appropriate conditions. In this case, when 
n = 3 the equation designs the physical KGD model developed by Geertsma and de Klerk [9 
and Khristianovich and Zheltov m- It represents the influence of the pressure exerted by a 
viscous fluid on a fracture in an elastic medium subject only to plane strain. This equation is 
derived from the conservation of mass for the fluid inside the fracture, the Poiseuille law and an 
appropriate pressure law (see [T5] section 3] and [H] for further details). In [T5] , weak solutions 
are constructed by passing to the limit in a regularized problem. The necessary compactness 
estimates are obtained from appropriate energy estimates. 
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The equation under consideration 


( 2 ) 


u t + d x (u n d x I(u)) = 0 

is a nonlocal degenerate parabolic equation of order a + 2. 

When a = 2, this equation coincides with the thin film equation (TFE for short) 

u t + d x (u n dl xx u) = 0. (3) 

This is a fourth order nonlinear degenerate parabolic equation originally studied by Bernis and 
Friedman [3] . This equation arises in many applications like spreading of a liquid film over a solid 
surface (n = 3) and Hele-Shaw flows (n = 1) (see [121UU [T2J [TS] El G3 [2]). TFE is derived also 
from a conservation of mass, the Poiseuille law (derived from a lubrication approximation of the 
Navier-Stokes equations for thin film viscous flows) and various pressure laws. The parameter 
n £ (0, 3] models various boundary conditions at the liquid-solid interface. The case n > 3 is 
mainly of mathematical interest m- In [3] weak solutions u are exhibited in a bounded interval 
under appropriate boundary conditions. In addition, they proved that u is nonnegative if uo 
is also so, and that the support of the solution u(t ,.) increases with t if uq is nonnegative and 
n> 4. 

For a = 0, the porous medium equation (PME for short) is recovered 

u t - d x {u n d x u) = 0. (4) 

This is a nonlinear degenerate parabolic equation. The simple PME model describes the modeling 
of the motion of a gas flow through a porous medium |20j. In this case, the PME is derived from 
mass balance, Darcy’s law which describes the dynamics of flows through porous media, and 
a state equation for the pressure m- PME also arises in heat transfer m and groundwater 
flow (19l and was originally proposed by Boussinesq. It took many years to prove that PME 
is well posed and the famous source type solutions were found by Zel’dovich, Kompanyeets 
and Barenblatt [20:. The questions of existence, uniqueness, stability, smoothness of solutions 
together with dynamical properties and asymptotic behavior are well represented in [20| where 
two main problems are studied. First, the domain space is K d and the initial condition uq has a 
compact support so the solution u(t, x ) vanishes for all positive times t > 0 outside a compact set 
that changes with time. Secondly, if the initial data has a hole in the support then the solution 
has a possibly smaller hole for t > 0. 

Note that TFE can be seen as a fourth order version of the classical PME [13]. Furthermore, 
both equations are parabolic in divergence form. In both cases, there are compactly supported 
source type solutions (n > 1 for PME (201 and 0 < n < 3 for TFE [4]) [3]. The most famous 
common properties are finite speed of propagation and the waiting time phenomenon. Similar 
properties are expected in our case. Self-similar solutions are constructed in El but other 
properties are still not proved. One striking difference between TFE and PME is the lack of a 
maximum principle for TFE 0. 

The case a £ (—2,0) corresponds to the fractional porous medium equation studied in S7]. 
Explicit self-similar solutions are exhibited and, under appropriate conditions, weak solutions 
are constructed. 

In this paper, we will generalize the result of [TSJ to the cases 0 < a < 1 and 1 < a < 2. We 
prove a result of existence with the same approach as that in the case a = 1 but by modifying 
the compactness results. Consequently all cases a £ [0,2] are now covered. 

In the case a > 1 we get the local uniform convergence of approximate solutions due to the 
following embedding in dimension 1 

H%(n) C'°’^(D). 
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This convergence allows one to pass to the limit in the nonlinear term and then allows us to 
construct nonnegative solutions for nonnegative initial data merely in H* (SI). 

In the case a < 1 because of the following embedding 

H%(fl) ^ L p (n) for all p < 

1 — a 

we can get a compactness result in L p (fl) only for p < and not for all p < oo as in the case 
a = 1. Neverthless, we recover a compactness result for the term I(u) which allows us to pass 
to the limit and conclude. 

In both cases, we prove that the solution is strictly positive under a condition on n. 


Integral inequalities 

Assume that SI = M, if u is a solution of m then it satisfies the energy inequality 
r rT r r 


— / u(t)I(u(t))dx + 2 / / u n d x I(u) 2 dxdt < — / uoI(uo)dx. 


J n J o Jci Jn 

Observe that — J ul(u) is the homogeneous H% norm. Let G be a nonnegative function such 
that G"(s) = 4r- Then the positive solution satisfies 


G(u(t))dx 


d x ud x I(u)dxdt ^ 


f G(uo)dx. 
J n 


Note that — / d x ud x I(u) is the homogeneous H£ +1 norm (it is in fact a Neumann-Sobolev space, 
see below). We see that the energy inequality controls the L°°(0, T;H%(£1)) norm of the solution. 
For the function G mentioned above, we can take 


G(s) 



(5) 


so that G is a nonnegative convex function satisfying G(l) = G'(l) = 0, G(s) = oo for all s < 0 
and for s > 0, we have 


s In s — s + 1 

when n = 1 

s 2 ~ n , a , 1 

(2— n)(n— 1) ' n— 1 ' 2 —n 

when 1 < n < 

h^ + s -l 

when n = 2 

1 1 | S 1 

when n > 2. 

(n—2)(n— 1 ) s n ~ 2 ' n —1 n— 2 


Main results 

In this work, we prove three main results. We first prove the existence of nonnegative weak 
solutions for the problem with 0 < a < 1 for nonnegative initial data with apropriate conditions. 
Secondly, for a > 1, we construct nonnegative solutions for nonnegative initial data in H~(Q). 
Finally, we prove the strict positivity of solutions for large n's. 

Theorem 1.1 (Existence of solutions for 0 < a ^ 1). Let n ^ 1 and a £ (0,1]. For any 
nonnegative initial condition uq £ H 2 (Jl) such that 


/ G(uo)dx < 00 
Jn 


( 6 ) 
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where G is a nonnegative function such that G"{s ) = -4tj there exists a nonnegative function 
u £ L°°(0, T; (17)) n L 2 (0, T; ff| +1 (17)) 
which satisfies on Q = (0, T) x 

// udtspdtdx — // nu n ~ 1 d x ul(u)d x ipdxdt — // u n I{u)d xx ipdxdt = — / uo<£>(0, .)dx (7) 

JJq J./q jJq in 

/or all ip £ 27([0,T) x 17) satisfying d x tp = 0 on (0,T) x dfl. 

Furthermore u satisfies for almost every t £ (0,T) 


u(t,x)dx= / uo(x)dx 
i dn 


( 8 ) 


and 




/ / g 2 dxds < ||uo||^a 

o do 

, — 2 


where the function g G L 2 (Q) satisfies g = d x (u 2 /(w)) — 2 d x ul(u) in V'(Q), and 

[ G(u(t,x))dx + [ ||n|| 2 « +1 ds < [ G(u 0 )dx. 

Jn Jo h n ( fi ) dn 


(9) 


( 10 ) 


Remark 1.1. The weak formulation o comes after two integrations by parts of the equation @. 
We recall that the function G : R + —> R + is given by ©■ Note that the space H^( 17) is defined 
via spectral decomposition of — (—A)a (see below). 

Theorem 1.2 (Existence of solutions for 1 < a < 2). Let n ^ 1 and a > 1. Tor any nonnegative 
initial condition uq £ (17) , 7/iere exists a nonnegative function 




(Q) 


such that 


d x I(u) £ L? oc (Q+) 


( 11 ) 


and that satisfies 



udtipdtdx + 



u n d x I{u)d x ipdxdt 


UQip(0, .)dx 


( 12 ) 


where Q+ = {n > OjnQ, for all ip £ 2?([0, T) x 17) satisfying d x ip = 0 on (0, T) x917. Furthermore, 
u satisfies conservation of mass. 

Theorem 1.3 (Strictly positive solutions). Assume 0 < a < 2 and n > 2 + TTiere exists 
a set P C (0,T) such that | (0,T) \ P |= 0 and the solution u constructed as in Theorem, II.II 
satisfies u(t ,.) £ C 0 ’^(17) /or allt £ P and for all /? < min{l, ^1} and u(t,.) is strictly positive 
in 17. Furthermore, u is a solution of 


ut + d x J = 0 in D'(i 7) 


where 


J(t ,.) = u n d x I{u) £ L 1 (17) /or all t £ P. 
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Organization of the paper 

The paper is organized as follows: in Section 2, we define the nonlocal operator I by using the 
spectral decomposition of the Laplacian and we write an integral representation for it. Then 
we prove two important Propositions used in the proofs. In Section 3, we study a regularized 
problem before proving our Theorems in Section 4. 


Notation 


In this work, we denote II = (0,1) and Q = (0 ,T) x SI. The space if^(SI) is the functional space 
defined in nu Section 3.1] by 

( OO CX) ^ 

H s N (Sl) = s u = ^ Ck<fik', 4(1 + X s k ) < +oo > 


fc=0 


fc=0 


where {X k , tp k }k>o are the eigenvalues and corresponding eigenvectors of the Laplacian operator 
in SI with Neumann boundary conditions on 9SI with the norm 


IMI/r« (n) — + ^fc)> 


equivalently to 


Mla^n) 


k =0 


udx] + l |w||2 ff ; ( n) 


where the homogeneous norm is given by 




= Y, 4 K- 


k =0 


Note that ff^-(SI) = H S (Q) for all 0 < s < § (see |T|) with equivalent norms. Indeed, 


L N 

IMI^(n) = IMlI’pi) + IMI ff>) 

and since we are in dimension 1 we have for these values of s 

Ml 


H e N ( n) 


Note also that we have 


f Q udx < C(SI)||u|| 2 (Holder inequality), 


Ml! ^ (7(11)11 (—A) 2 u||| < c||u|| 2 a (fractional Poincare’s inequality). 

H n (Q) 


Finally, as usual s + = max{0, s}. 


(13) 


(14) 


(15) 


2 Preliminaries 

2.1 Operator I 

Spectral definition. We define the operator I by 

OO OO 

I : c k tp k — > — c kX k (p k which maps H^(SI) onto L 2 (SI) 

fc=0 fc=0 
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where {A&, tfk}k >o are the eigenvalues and corresponding eigenvectors of the Laplacian operator 
in SI with Neumann boundary conditions on dfl: 


{ -A ip k = \ k ipk in SI, 
d v (fik = 0 on 1, 

In ( Pk dx = 1 - 


Integral representation. The operator / can also be represented as a singular integral oper¬ 
ator. We will prove the following. 

Proposition 2.1. Consider a smooth function u : SI — > R. Then for all x € fl, 

I{u)(x)= / (u(y) — u(x))K(x, y)dy 
J n 

where K(x,y) is defined as follows. For all x,y € SI 


K(x,y) =c a ^2 


fee i 


1 


1 


x-y-2k\ 1+a | x + y - 2k | 1+c 


where c a is a constant depending only on a. 

Proof. Let’s replace SI by (—1,1) and u by its even extension to (—1,1). Then let’s extend u 
periodically to R and let u be this extension. For x £ SI, 


I(u)(x) = -(-A) 2 u(x) = c a [ ( u(y ) - u(x)) 

Jr 

p l+2fc 

= c «m/ (u{y)-u( x )) i 


y — x | 1+Q 


dy 


-l+2fc 


y ~ x | 1+a 


= c a (u{y) - zt(z)) V 

-'O 


— 7? |!+ a 


dy 


— y — 2k | 1+a | x + y — 2k | 1+0 


because u is 2-periodic 

because u is even. 


□ 

Now we can easily conclude the following Corollary. 

Corollary 2.1. Consider two smooth functions u, (p : S7 — > R. Then 

/ I(u)(x)ip(x)dx = / u(x)I(tp)(x)dx (16) 

J f2 J 

2.2 Important identities 

As [151 Section 3], the semi-norms ||.||^ , ||-|| ff “ (a) , |M|^f+i and ||.||^«+i (fJ) are related to 

the operator I by important and very useful equalities. 

Proposition 2.2. 1. For all u £ H%(£!), we have —(I(u),u) = llull 2 * 

h 2 (n) 
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2. For all u £ Cl), we have ||w|| 2 a = J Q I(u) 2 dx. 

H n (0) 


3. For all u £ ff^ +1 (fl), we have ||u|| 2 a+1 = — f n I(u) x u x dx. 

Hn (O) 1 

4■ For all u £ i?^ + 1 (fi), we have ||u|| 2 a+i = J Q I(u) x dx. 

h n (fi) 

Proof. Note that if u £ H% (Cl) then I(u) £ H~% (f 1) and 

OO 

< I (u),v) H - 9(n)> H9 (n) = -'52 c k AJ4 

k =0 

where v = XlfcLo dkVk G H% (fi) and u = c fc<Pfc> so 


- uUu) = £ cW = |M^, 

k =0 

The second equality is actually very easy to prove since I(u) = — YlkLo CkP£ ■ Indeed, 


i{u) 2 dx = y 


fc=o 


« = 11 - 11 ^( 0 )' 


In order to prove the other equalities, we note that (d x (p k )k form an orthogonal basis of 


L 2 (fl). We write 


u x = ^ ~2c k d x (fk in L 2 (fl). 
k =o 

OO 

and d x I(u) = - ^ c k X k <9^ in L 2 (f^) 


fc=i 


so 


n OO n OO n 

- / I{u) x u x dx = ^2 c l X k / dxpidx = ^2 c l X k / fk{-d xx ipk)d ; 
fc=0 fc=o 

OO 


= H C k X k / x kVkdx = J2 c t X k + = 


k—0 


k—0 


H£ (fi) 


For the last equality, 


OO « OO 

/ I(u)ldx = J2 c lK / d x ip 2 k dx = y] c feA“ +1 = IMI^=+i ■ 

fe=l ■' n fe=0 N ' 


□ 


2.3 The problem —I(u) = g 

We consider the following problem 

J For a given 3 £ L 2 (f2), find u £ i7^(f2)such that 
1 -I(u)=g. 


(17) 


Since f Q I(u)dx = 0 for all u £ fl), we must assume that f Q g(x)dx = 0 otherwise (1171) 
has no solution. 
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Proposition 2.3. For all g £ L 2 (Fl) such that J^gdx = 0, there exists a unique function 
u £ such that 


—I{u)=g in L 2 (Q) and / udx = 0. 

Jn 


Furthermore if g £ H 1 ^), then u £ iL)0 +1 (fi). 

Proof. Let g £ L 2 {Pi). For g = YJJk=i dkPk with J2kLi d1 < °°> we consider 

, ^ d k r 

u = I (g) = 2 _^ S FFjf(Fl) and verify / udx = 0. 

fc=l ^n 

Since (ip k ) k form an orthogonal basis of L 2 (fl), the solution is the unique satisfying j n udx = 0. 
It is clear that every further regularity on g will imply a further regularity on u shifted by an 

a. □ 

We thus conclude the following Corollary which will be used to prove the existence of solutions 
for the stationary problem. 


Corollary 2.2. For all g £ L 2 (Ft), there exists a unique function v £ F[ff(fl) such that 


— I(v) + / vdx = g. 
J n 


(18) 


Furthermore if g £ H 1 (f2), then u £ iL^ +1 (f2) and the map g —► u is bijective. 

Proof. Let m = j n gdx and g' = g — m. Then g' £ L 2 (fl)(since fi is bounded) and g'dx = 0. 
From Proposition 12.31 there exists a function u £ F[^(fl) such that 

—I(u) = g' and / udx = 0 . 

Jn 

Let v = u + m. Then vdx = m and 

—I{v) = —I(u) =g' = g — m = g— / vdx. 

Jn 

For the uniqueness, consider two solutions v\ and V 2 then 

/ Vidx = / V 2 dx = / qdx 

Jn Jn Jn 

and w = V\—V 2 satisfies — I(w ) = 0. Hence, w = 0 from the uniqueness given by Proposition ^. 31 

□ 


3 Regularized problem 


We consider the following regularized problem 


I d t u + d x (f e (u)d x I(u)) = 0 
< d x u = 0, f € (u)d x I(u) = 0 
[it (a;, 0) = u 0 ( x) 


for x £ fl, t > 0 , 
for x £ t > 0 , 

for x £ SI, 


(19) 


8 


where / £ (s) = s" + e, e > 0 and 0 < a < 2. 

To prove Theorem ll.ll and ll.2l we need to prove the existence of a solution for the regularized 
problem. Let us pass with the following stationary problem 


For t > 0,g £ H 2 (fi), find u € s.t. 


u + Td x (f e {u)d x I(u)) = 5 
d x u = 0 ,d x I(u) = 0 


in Q, 
on (9S2. 


( 20 ) 


Once we get a solution for (Em> , we can prove the existence of a solution for 


3.1 Stationary problem 

Proposition 3.1. For all g £ there exists u £ such that for all ip £ H 1 ^) 

we have 



— T 


f e (u)d x I(u)d x ipdx 


/ gydx. 
J n 


Furthermore, u verifies 


( 21 ) 


u{x)dx = / g(x)dx 

i Jn 


( 22 ) 


and 


i|| 2 f +2 t f f e (u)d x I(u) 2 dx ^ || S || 
h T (n) J n 


n) 


(23) 


If f n G e (g)dx < oo where G £ is a nonnegative function such that G"(s) = Yfsj’ then 


[ G e (u)dx + t||? 

Jn 


I 2 c 

FT ' 

■tl ^ 


(O) 




Ge(g)dx 


(24) 


Remark 3.1. Note that we can consider G £ (s) = /* G'f(r)drdt , so G £ is a non- negative convex 

function for all e > 0 satisfying G £ (l) = G'(1) = 0. 

Proof. Thanks to Corollary [221 we can recover all test functions from H 1 (n) by considering 


<£> = —I(v) + / ud:r 

Jn 


for some function v £ Hf[ +1 (n). So equation (T2Tt becomes 


/ ul(v)dx +( / udx I I / ud:r I + r / f e (u)d x I(u)d x I(v)dx 
Jn \Jn ) \Jn ) Jn 

= — ql(v)dx + ( / qdx ) I / vdx 1 . 

Jn \Jn J \Jn J 


(25) 


Now, we consider the nonlinear operator A defined by 


A(u)(v) = — ( ul{v)dx +( [ udx I f f vdxj +t f f e (u)d x I(u)d x I(v)dx for u, v £ U^ +1 (n). 
Jn \Jn J \Jn J Jn 
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We prove that this is a continuous, coercive and pseudo-monotone operator. Note that the 
functional T g defined by 


Tg{v) = — [ gl(v)dx + ( f gdx] f f vdxj for v £ H^ +1 ( fi). 

Jn \Jn ) \Jn ) 

is a linear form on H c ^ +1 ( fi). So our problem reduces to the following 

{ Let V = H« + \n). 

A : V —>■ V* coercive, continuous and pseudo-monotone. . 

T, e V. <26) 

Find u £ H^ +1 ( fl) such that A(u ) = T g in V*. 

The theory of pseudo-monotone operators m implies the existence of a solution for m so there 
exists u £ such that 

A(u)(v) = T g {v) for all v £ H^ +1 (n). 

It remains to prove that A is a continuous, coercive and pseudo-monotone operator on H^ +1 ( fl). 
The reader can find the proof in H3 Appendix A] for V = H^( f2) but this proof can be easily 
adapted for our case V = H^ +1 ( fi). 

By using Corollarv l2.2l we deduce that u satisfies (1211) for all ip £ H 1 ( fl). 

For the properties of u. first by taking <p = 1 as a test function in m we obtain mass conser¬ 
vation (1221) . Secondly, take v = 1 1 — Jq udx in (ESI), by using Proposition 12.21 we have 


i|| 2 f -Hr / f e (u)d x I(u) 2 = — f gl(u)dx 
H^(n) J n J Q 


€ 


'H* 


1 


< ollffll 


(1) 


+ dM 


H^( fi) 


which ESI (Note that the high regularity of g is solely used in this inequality, otherwise g £ L 2 (fi) 
is sufficient to prove the existence above). Finally, note that G' e is smooth with G' e and G" are 
bounded, and Q is bounded so we can take p = G' t (u) £ as a test function in ED' 


[ uG' t {u)dx — t j f e (u)d x I(u)d x uG"(u)dx = f gG' e (u)dx. 

J Cl J Cl J Cl 

So by using Proposition 12.21 and the fact that G"(s) = we get 

t || u || 2 f+1 = f G' e (u)(g — u)dx < f (G e (g) - G e {u))dx 

h n (°) Jn Jn 

because G e is convex and we deduce El- 


□ 


3.2 Implicit Euler scheme 

We construct a piecewise constant function 

u T (t, x) = u k {x) for t £ \kr, (k + 1 )r), fc £ {0,..., N — 1} 
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where r = and (M fc )fcg{o....,Ar-i} is such that 

U fc+1 + T d x (f e (u k+1 )d x I(u k+1 )) = u k . 

The existence of the u k follows from Proposition 13.11 bv induction on k with u° = u 0 . We deduce 
the following 

Corollary 3.1. For any N > 0 and Uq £ (12), there exists a function u T £ L°°(0,T; (12)) 

such that 


1. t —>■ u T (t , x) is constant on [hr, (k + 1)t), k £ {0 ,N — 1} and r = jj-. 

2. u T = uq on [0 ,t) x 12. 

3. For all t £ (0, T), 


/ u T [t 1 x)dx= / uo(x)dx. 

Jn J n 


(27) 


l For all ip £ (7,1(0,T; H 1 (fl)), 

u T — S t u t 

J Qt,T 


T JJQ 

where S T u T (t,x) = u T (t — t,x) and Q t ,t = (t, T) x 12. 
5. For all t £ (0, T), 


ipdxdt = // f e (u T )d x I(u T )d x ipdxdt (28) 

J t 


IKCMIllf 


H 2 (n) J o J n 


f t {u T )d x I{u T ) 2 dxdt < ||m 0 || 2 f 


h 2 (n) 


6. If f n G e (uo)dx < oo, then for all t £ (0, T) 


(29) 


[ G e (u T (t,x))dx+ [ |K(s,.)|| 2 f+I ds < [ G e (u 0 )dx. (30) 

Jn Jo h n (°) Jn 

3.3 Existence of solution for the regularized problem 

Now we are able to prove the existence of a solution for the regularized problem. 

Proposition 3.2. Let 0 < a < 2. For all uq £ (11) and for all T > 0, there exists a function 

u e such that 


satisfying 


u £ £ L“(0,T;tft(12))nL 2 (0,T;i7“ + 1 (12)) 


u e dtpdxdt 


Q 



f e (u e )d x I(u e )d x ipdxdt 


zio</?( 0 ,.)dx 


for all ip £ C 1 ( 0 ,T; i7 1 (n)) with support in [ 0 ,T) x 12 . 

The function u e satisfies for almost every t £ (0 , T) 

/ u e (t,x)dx = / uo(x)dx 

Jo, Jo, 


(31) 


( 32 ) 
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and 


||i/(t,.)|| 2 * +2 f j f e (u e )d x I(u e ) 2 dxdt < |K|| 2 « . (33) 

Finally, if f n G e (uo)dx < oo then for almost every t £ (0 ,T), 

[ G e (u e (t, x))dx + [ ||u £ (s, .)|| 2 f+i ds < [ G e (u 0 )dx. (34) 

J n Jo h n ( n ) Jn 

Proof. We consider the sequence ( u T ) constructed in Corollary 13.II and let r —> 0. Bound (E51) 
and (B3) implies that ( u T ) is bounded in L°°(0,T; (12)) and ( d x I(u T )) is bounded in L 2 (Q). 

Case 0 < a ^ 1. Note that 

MT ~ ^ = d x (fe(u T )d x I(u T )). 

T 

Since n > 1, the function f e is Lipschitz and so ( f t {u T )) is bounded in L°°(0, T thus by 
the Sobolev embedding theorem, we deduce that ( f e (u T )) is bounded in L°°(0, T; L p (12)) for all 
p < jz^- We know that (■ d x I{u T )) is bounded in L 2 (0, T; L 2 (12)) so f e (u T )d x I (u T ) is bounded in 
A 2 (t,T; 1/(12)) where \ = \ + y We deduce that 

d x (fe(u T )d x I{u T )) is bounded in L 2 (t, T; W _ 1 ’ t ’(12)) 

Since a ^ 1, we have the following embedding 

H*(Q) L p (12) ->■ W _:M (12) 

for all p < and for all2 > 2 (because 12 is bounded and we have a Sobolev space of negative 
regularity). Aubin’s lemma implies that ( u T ) is relatively compact in C°(0,T; 1/(12)) for all 
p < Note that ( d x I(u T )) is bounded in L 2 (12) and (it r ) is bounded in L°°(0, T; 2/(12)) 

(because 1 < jz^)- Hence, ( u T ) is bounded in L 2 (0,T; i?^ + 1 (H)). Since 

tf“ + 1 (H) ^ 22 | + 1 ( 12 ) W _ 1 ,i ( 12 ), 

we deduce that ( u T ) is relatively compact in 2/(0, T; 22^ + 1 (12)). So we can extract a subsequence, 
also denoted ( u T ), such that when r tends to zero we have 

• u T —> u e £ L°°(0,T; (fl)) almost everywhere in Q, 

• u T —> u e in 2 /( 0 , T; 22 ^ + 1 ( 12 )) strongly, 

• d x I(u T ) — 1 d x I(u e ) in L 2 (Q) weakly. 
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Now let us pass to the limit in (051) . We have 


u T - S t u t 1 

- <pdxdt = — 

Qt,T t t 


u T (t, x)dtdx — / / u T (t, x)ip(t, x)dtdx 


o Jn 

T-t r 


io Jn 


u T (t, x)if(t + t , x)dtdx 


/ o Jn 

Jo Jn T 

T r 


1 


1 


u T {t , x)<p(t, x)dtdx H— 


u T (t, x)tp(t + t, x)dtdx 


1 0 Jf2 


—>— // u e d t pdxdt— / u £ (0,a;)<p(0,x)da; + 0. 
T_! '° JJq Jn 

For the nonlinear term, we integrate by parts 


f e (u T )d x I(u T )d x <p = - / I f e (u T )I(u T )d 2 xx <p- // n{u T ) n 1 d x u T I(u T )d x (p. 

'q JJq 


u T —> u e in L 2 { 0, T ; for all s < 1 + a. 


I{u T ) I(u e ) in L 2 (0,T;H s \n)) for all s' < 1 


We have 


So 


and 


d x u T —> d x u e in L 2 (0 ,T]H S (f Vj) for all s" < a. 
So we deduce the following convergences 

I(u T ) ->■ I{u e ) in L 2 (0,T;L 9 (tt)) for all q < oo. 
u x —>■ u% in L 2 (0,T;L p (fI)) for all p < 


1 — a 


Furthermore, since u T u e in C'°(0, T; L p (fl)) for all p < and f e is lipschitz then 

2 


f t {u T ) f e {u e ) in C u (0, T; L p (fl)) for all p < 


1 — a 


For the term ( u T ) n - 1 , if n ^ 2 then the function s —> s n 1 is lipschitz and 


(u T ) n - 1 (u e ) n in C u (0,T-,L p (Cl)) for all p < — 


a 


If n < 2 then —^-r ^ 1 and 

n— 1 


(35) 


(u T ) n ~ l (it 6 )"" 1 in C u (0 ,T;L— (ft)) for all p < 


1 — a 
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Thus we can pass to the limit in (1351) and reverse the integration by parts to obtain 



f e (u T )d x I{u T )d x ip 


fe{u e )I(u e )d xx ip - // n(u e ) n 1 d x u e I{u e )d x ip 


Q 


j j f e (u e )d x I(u e )d x ip. 


For the properties of u e , first since u T u e in L°°(0, T; L 1 (fl)) mass conservation equation (1551) 
follows from m- 

Secondly, we note that ( u T ) is bounded in L°° (Q,T; H% (£1)) so ( u T ) weakly converges to u e 
in H * (SI) and 


'h 2 (n) 


< liminf ||u r || . a 

“ t^O 11 n H 2 (Q) 


Note that estimate (l29l) implies that y/ f e (u T )d x I ( u T ) is bounded in T 2 (0, T; L 2 (£7)) thus it weakly 
converges in L 2 (0,T; L 2 (Vl)) and the lower semicontinuity permits us to conclude (1551) . 

Finally, to derive (1551) we note that G e (u T ) —> G e (u e ) almost everywhere and Fatou’s lemma 
implies for almost every t £ (0, T) 


GJu e (t,x))dx < liminf [ G e (u T (t,x))dx. 
! J n 


Furthermore, ( u T ) is relatively compact in L 2 (0, T; Hj^ +1 (n)) thus 


/ ||M e (s)||2 f+1 ds = lim f ||u r (i 
Jo h n T_>u d 0 


Hi 


_%+ids. 


Hence (1501) implies (1551) . 
Case 1 < a < 2. Note that 


uT StuT = d x (f e (u T )d x i(u T )). 

T 

We have ( u T ) is bounded in L°°(0, T ; (S2)) so by the Sobolev embedding theorem, we deduce 

that ( u T ) is bounded in L°°(0, T; C 0 ’^ (fi)). Thus (/ e (w T )) is bounded in L°°(0, T; L°°(fl)). We 
know that (d x I(u T )) is bounded in T 2 (0, T; L 2 (S2)) so ( f e {u T )d x I(u T )) is bounded in L 2 (r, T; L 2 (fl)) 
We deduce that 


d x {f e (u T )d x I(u T )) is bounded in L 2 (t,T\W 1 ,2 (S2)). 

Since a > 1 we have the following embedding 

H%(n) ^ C°’^(H) -> w~ h2 {n). 

Aubin’s lemma implies that the sequence ( u T ) is relatively compact in C°(0, T; C 0, ~~ (f2)). Since 
(■ d x I(u T )) is bounded in L 2 (H) and ( u T ) is bounded in L°°(0,T; L 1 )^)) then, ( u T ) is bounded in 
L 2 (0,T; Using the following embedding 
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we deduce that ( u T ) is relatively compact in L 2 (0, T; So for a subsequence we have 

• u T —> u e locally uniformly, 

• d x I(u T ) — i d x I(u e ) in L 2 (<3)-weakly, 

• u T —> u e in L 2 (0,T; H^ +1 (fi)) strongly. 

Let us pass to the limit in (E51) . As in the first case 


-Sr 


JQt,T 

For the nonlinear term, since 


-pdxdt 


r—>0 


u e d t tpdxdt — / zt e (0, x)tp(0, x)dx. 


u T —> u e locally uniformly, 


Then 

f e (u T )d x ip f e (u e )d x tp in L 2 (0,T; L 2 (H)) - strongly. 

Furthermore 

d x I(u T ) —- d x I{u e ) in L 2 (0, T; L 2 (Q)) — weakly. 

Hence 



f e (u T )d x I(u T )d x ipdxdt -)• 



f e (u e )d x I(u e )d x tpdxdt 


and the proof is complete. 

For the properties of u e , the proofs of estimates ©, (El and are the same as in the 
first case. □ 


4 Proofs of main results 

4.1 Proof of Theorem 11.11 

Consider the sequence (w e ) such that u e £ L°°(0, T; (H)) fl L 2 (0, T; solution of 

m- Our goal is to pass to the limit e —> 0. 

Note that (1351) and (15^1) imply that ( u e ) is bounded in L°°(0, T; (fi)). Since f e is Lipschitz 
then f e (u e ) is bounded in L°°(0, T; (fi)). So by using the Sobolev embedding theorem, we 
deduce that / e (u e ) is bounded in L°°(0, T; L P (Q,)) for all p < 

Furthermore, fITTTTT) also implies that f e (u e )?d x I(u e ) is bounded in L 2 (0, T; L 2 (fi)). Thus 

f e (u e )d x I(u e )is bounded in L 2 (0,T; L r (£l)) where - = - + —. 

r 2 2p 

Hence 


d t u e = —d x (f e (u e )d x I(u e )) is bounded in L 2 (0,T; W 1,T ’(H)). 


Since 


H*(Q) (12) -► 
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Aubin’s lemma implies that ( u e ) is relatively compact in C°(0, T; L p (fl)) for all p < So we 
can extract a subsequence such that 

• u e ^um C°(0,T;LP(n)) for all p < j^. 

• u e —> u £ L°°{ 0, T; (Q)) almost everywhere in Q. 

Let us pass to the limit in (15T1) . Let tp € D([0,T) x fi) satisfying d x p = 0 on (0,T) x dfl. 
Since u e — > u in C°(0, T; L 1 (n)), we have 


Q 


u e dt<pdxdt —> / / udtpdxdt. 


Q 


Remark that (1551) implies that 


e JJ ( d x I{u e )) 2 < c. 

The Cauchy-Schwarz inequality implies 

e JJ d x I(u e )d x ipdxdt < c(p)v^(v / e||3 :E /(w e )||2) ->• 0. 


Estimate (1551) also gives that (u e )+d x I(u e ) is bounded in L 2 (0,T; L 2 (S1)). For the term (u e )%, 
we consider two cases, if n ^ 2 then the function s —> s * is Lipscliitz and ((u e ) * ) is bounded in 
L°°(0,T;L p (fl)) for all p < j ry- We deduce that ((u £ )” d x I(u e )) is bounded in L 2 (0, T; L m (fl)) 

where ~ = ^ If 7Z < 2 then ((it e )^) is bounded in L°°(0,T; (f2)) for all p < (in this 

case pp ^ 1). We deduce that ((u e )™d x I(u e )) is bounded in L 2 (0,T; L m (Q)) where ^ | 

hence 


h e := (it e )™ d x I(u e ) — 1 h in L 2 (0,T; L m (Tl)) weakly. 
Passing to the limit we obtain 

/ / udt<pdxdt + / / hd x (pdxdt = — / uo<p(0,x)dx. 

J Jo J Jo J n 


'Q 

It remains to show that 


/i = u+d x I(u) 


in the following sense 

// hipdxdt = — // nu r ^ l d x ul(u)ipdxdt — // u+I(u)d x ipdxdt 

JJq JJq JJq 

for all test functions <p such that tp = 0 on (0, T) x <9fl, that is 

/i = d x (u+I(u)) — nu r J t T 1 d x uI(u) in D'(Q). 

Note that G e is decreasing with respect to e, so 


/ G e (uo)dx < / G(uo)dx < c. 

J n 


(36) 
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Thus estimate O implies that ( u e ) is bounded in L 2 (0, T; H£ +1 (Q)). Recall that ( d t u e ) is 
bounded in L 2 (0 ,T;W~ 1,l (fl)). Aubin’s lemma implies that 

(■ u e ) is relatively compact in T 2 (0,T; for all s < ^ + 1. 


Hence 


( d x u e ) is relatively compact in L 2 (0,T- 1 H s (H)) for all s' < ^ 


and 

(/(u e )) is relatively compact in L 2 (0,T; (fl)) for all s" < 1 — —. 

Thus we can extract a subsequence such that 

u e —>■ u in C°(0,T;L p (Q)) for all p < - -, 

1 — a 

I{u e ) — > I(u) in L 2 (0,T; L q (Q)) for all q < oo, 

d x u e —> d x u in L 2 (0,T; L P (Q)) for all p < —-—. 

1 — a 


We write 



h e ipdxdt 


( u e ) r i_d x I(u e )(pdxdt 



d x u e I{u e )ipdxdt 


(u e )1I(u e )d x <pdxdt. 


Using these convergences and the fact that I(u e ) converges in L 2 (0, T; L q (fl)) for all q < oo we 
can pass to the limit and obtain (1331) . Note that for the terms ( u e ) n and (u e ) ra_1 we consider 

two cases n ^ 2 and n < 2 and we proceed as above. In the first case the functions s —> s n 

and s — > s n ~ 1 are Lipschitz and then ( u e ) n — > u n and (u e ) ra_1 u n ~ l in C°(0, T; L p (fl)) for 

all p < If n < 2 then £ > 1 and ^ 1 and ( u £ ) n — > u n in C°(0, T; L% (fi)) and 

{u e ) n ~ l —> u"" 1 in U°(0 ,T;L^t(U)) for all p < 

For the properties of u, passing to the limit in (1321) implies mass conservation equation (0. 
Since ( u e ) is bounded in L 2 (0,T; then u e —*■ u and 


1 l 2 (o,t-,h\ 


i+i <liminf| 

r m - e—>0 1 


m’ 


Note that 


G e (u e ) G{u ) almost everywhere and G e (uo) < G(uo). 

Then by Fatou’s lemma estimate m follows from (1331) . 

Remark that estimate (1331) implies that g e = (u e )ld x I(u e ) weakly converges in L 2 (Q) to a 
function g and the lower semi-continuity of the norm implies 0. It remains to prove that 

n fl XL _ 1 

g = d x (u+I(u)) - -u| d x ul(u) in D'(Q). (37) 
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We have 



g e <pdxdt 


JJ ( u e )+d x I(u e )ipdxdt 
— JJ 1 d x u e I(u e )tpdxdt 


(u e )+I ( u e )d x ipdxdt. 


Also, using the convergences above and the fact that I(u e ) —> I(u) in L 2 (0, T; L 9 (ft)) for all 
q < oo, we can pass to the limit and obtain (EH). Note also that for the terms (it e )2 1 and (u e )i 
we consider two cases n ^ 4 and n < 4 and we proceed as above. 

It remains to prove that u is a nonnegative function. Note that estimate (1341) implies that 
for all t e (0 ,T) 


Since 


G e (u e (t : x))dx < / G e (uo(t,x))dx. 


we conclude that 


G e (uo(t, x))dx < / G(uo(t, x))dx < oo, 


Note that for all 5 > 0, 


limsup / G e (it e (f, x))dx < oo. 
e —>-0 J £7 


(38) 


lim G e (—<5) = +oo. 

e-s-0 


Recall that u e {t,.) converges almost everywhere. So for 77 > 0, Egorov’s theorem implies the 
existence of a set A v c ft such that 

I ft \ A |< 77 and u e (t,.) —► u(f,.) uniformly in A v . 


Let S > 0. We consider 


C v> s = A v n {u(t ,.) < -26}. 

For every 77 , 6 > 0, there exists € 0 ( 77 , <5) such that if e < eo( 7 ?> <5) then rt e (f,.) < in C v j- 
This implies that C v j has measure zero. Indeed, if not then for e < £ 0 ( 77 , (5) we have 

G e (u e (t,x)) > G e (—6) —> + 00 . 

By Fatou’s lemma 

liminf / G t {u e {t,x))dx ^ / lim inf G e (u e (t, x))dx = +00 
€—>-0 j I €—^0 

which contradicts (1551) . 

Hence for all <5 > 0 and all 77 > 0, we have 

| {«(£, •) < -25} |<| G„,{ I + I ft \ A,, |< 77 . 

Thus, | {u(t ,.) < —26} |= 0 for all 6 > 0. We conclude that 

{«(£>•) < 0 } = (J 

k> 1 

has measure zero and so u(t : x) > 0 for almost every x £ Q and for all t > 0 . 


{«((,.)< ^ 4 } 
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4.2 Proof of Theorem 11.21 

We organize this proof in two stages. In the first stage we consider nonnegative uq € 
satisfying d6j and we prove the existence of solutions as in Theorem 11.11 In the second stage we 
use this information to prove the existence of solutions for nonnegative initial data which belongs 
to H% (f2). 


First stage Consider the sequence ( u e ) such that u e £ L°°(0,T; H? (fi)) n L 2 (0, T; H^ +1 (fl)) 
solution of di). Our goal is to pass to the limit e —> 0. Note that (ITTTT1) implies that ( u e ) is 
bounded in L°°(0,T\ H% (fl)). So by using the Sobolev embedding theorem, we deduce that (zz e ) 
is bounded in L°°(0,T; C' 0, ' S 2 i (f2)). Hence (/ £ (zz e )) is bounded in L°°(0, T; L°°(fl)). Furthermore 
(1551) gives that (/ £ (zz £ )5 d x I[u £ )) is bounded in L 2 (0,T; L 2 (H)). We deduce that 

(. f e (u e )d x I(u e )) is bounded in L 2 (0, T; L 2 (Q)). 


So 


d t u e = —d x (f e (u e )d x I(u e )) is bounded in L 2 (0,T; W 1,2 (H)). 

Since 

Aubin’s lemma implies that (zt £ ) is relatively compact in C°(0, T ; C 0, ^~ (f2)). So we can extract 
a subsequence such that 

u e —> u locally uniformly. 

Now let us pass to the limit in (EfTl) . Proceeding as in the case 0 < a < 1 we get the same results 
but it remains to prove the equation on h i.e. HMD- Since 

/ G e (uo)dx < / G(uo)dx < c, 

Jo, Jq 

estimate (151) implies that ( u e ) is bounded in L 2 (0, T; i?^ +1 (H)). Recall that ( d t u e ) is bounded 
in L 2 ( 0, T; fW 1,2 (f2)). Aubin’s lemma implies that 

( u e ) is relatively compact in i 2 (0, T; U S (H)) for all s < ^ + 1. 

Hence 

( d x u e ) is relatively compact in L 2 (0,T; H s (f2)) for all s' < — 

and since a < 2 then (m £ ) is relatively compact in L 2 (0, T; H a (Q)) and 
(/(z/)) is relatively compact in L 2 (0, T; L 2 {Vl)). 

Thus we can extract a subsequence such that 
z/ e —u locally uniformly, 

I{u e ) I{u) in T 2 (0,T;i 2 (H)), 

d x u e d x u in L 2 (0, T; locally uniformly with respect to x). 

Using an integration by parts for the equation of h e and using these convergences we can pass 
to the limit and obtain HMD- 

For the properties of u, the proofs are the same as in the case 0 < a ^ 1 but we use these 
convergences above to obtain the equation on g. 

We prove also that u is a nonnegative function as in the case 0 < a ^ 1. 
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Second stage Now we consider the case where u o ^ 0 belongs to H % (Q) without the additional 
condition ©. If we define 

u 0 s(x) = u 0 (x ) + S 

and denote us the nonnegative solution u constructed in the first stage for the initial data uos > 
which satisfies ©, then us satisfies 


us \< A , 


Q 


ugd x I(us) 2 dxdt < C, \ us(t,X2) — us(t,x\) \< k \ X2 — xi \ 2 , (39) 


with constants C , A, K independent of S and T. 

Proposition 4.1. There exists a constant M independent of 6 and T such that 

| u s (t 2 ,x) - us{ti,x) |< M | t 2 - h 1 2 <“+ 2 > (40) 

for all x € fl, ti and t 2 € (0 ,T). 

Proof. The proof is given in Appendix A. □ 

Taking a subsequence 

us — y u locally uniformly in Q, 

we will prove Theorem 11.21 Let <p G D([0,T) x Q) satisfying d x <p = 0 on (0,T) x dfi. We have 

// usd t (fdtdx + // usd x l(us)d x ipdxdt = — / (uq + 5)<p(0, .)dx. (41) 

JJq JJq Jn 

Since us —> u locally uniformly then 

/ / usdfipdtdx — > / / ud t ipdtdx and / (uo + <5)^(0, .)dx —> / mo</ 5 (0j 0^ as S —> 0. (42) 
JJo JJo Jn J n 


It remains to pass to the limit in the nonlinear term. We consider 

hs = Ugd x I{u s ). 

From (15(11) . ((us)%d x I(us)) is bounded in L 2 (Q) and since (us) is bounded in L°°(Q) so (hs) is 
bounded in L 2 (Q) and weakly converges to h in L 2 (Q). Our aim is to prove that 


h = 


u n d x I(u) in Q + := {it > 0} D Q 
0 elsewhere. 


For any r) > 0 we have 


O / ugd x I(us ) 2 > ( 7 :) 

J {u^rj}nQ 


' {u>r)}nQ 


d x I(us ) 2 


so ( d x I(us )) is bounded in L 2 ({u > 77 } D Q). Thus for all k € N, 

(d x I(us)) weakly converges in L 2 (Qk) 
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where Qk := {u ^ £} D Q. So, up to a subsequence, 

d x I(ug)) weakly converges to g in L 2 oc {Q + ) 

where <5+ = U Pk = {u > 0} fl Q. This implies that 
fcgN 

d x I(us) —*■ g in D'(Q + ). 


It remains to prove that 


g = d x I(u ) in D'(Q + ). 

Since ug —> u locally uniformly in Q then by using Corollary 12. II 

I {us) -t I(u) in V'{Q). 

So, d x I{ug) —> d x I{u) in V{Q). Now, let f £ V{Q + ) we have 

{d x I{ug),f)-D'(Q + )V{Q + ) = {d X I{ug),f)-D'(Q)T>(Q) {dxI{u),f)-D'(Q)-D(Q) 

where f is the extension by 0 of f to Q. So 

g = d x I{u) in V{Q+) and d x I{u) £ L 2 oc {Q + ). 

On the other hand, if 6 is sufficiently small, then 


f f ugd x I{u s )d x f |< c6?{ f f ugd x l{us) 2 )i < CS 

JJ{ u =o}nQ JJ 


(43) 


(44) 


Taking S —> 0 in (14T1) and using (l42ll , (14311 and (l44l) we deduce that m is satisfied. Finally since 
ug satisfies mass conservation and uniformly converges to u then u inherits the same property. 


4.3 Proof of Theorem 11.31 

Consider the sequence ( u e ) such that u e solution of (fT!)l) introduced in the proof of Theorem ll.il 
Recall that (28) implies that ( u e ) is bounded in L 2 {0,T-, H^ +1 {fl)). 

Case 0 < a < 1. We recall that (<9 t it e ) is bounded in L 2 {0, T; IT _1,i (f2)) . So Aubin’s lemma 
implies that ( u e ) converges in L 2 (0, T ; C°’ /3 (ri)) for all 0 < SL ^- We can thus find a subsequence, 
also denoted ( u e ), and a set P C (0, T) such that | (0, T) \ P | = 0 and for all t £ P, u e (t, .) 
converges strongly in C&{ fl). 

We note that for all t £ P, u is strictly positive. Indeed if there exists {to,Xo) £ P x 11 such 
that u(to, xq) = 0 then for any 0 < there exists a constant cp such that for all x £ 

u(t 0 ,x) < Cp | X - X 0 f . 


Thus 


G{u{to, x))dx > 


(cp\x- Xq I Py 


rdx. 
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Given n > 4, we can choose /? < such that /3(n — 2) > 1. We deduce 

J G(u(x, to))dx = oo 

which contradicts (1551) . 

We deduce that there exists 5 > 0 (depending on t) such that for e small enough 


u e (t ,.) ^ 6 in fl. 


Note that 


liminf / f e (u e ) | d x I{u e ) | 2 dx < oo for all t £ P. 
e ^° Jn 

Indeed, if we denote 

Ak = {t £ P; liminf / fAu e ) I d x I(u e ) I 2 dx ^ fcl 

Jn 

then using (1331) and Fatou’s lemma we have 

c ^ liminf [ [ f e (u e ) I d x I(u e ) I 2 dxdt 

Jo Jn 

^ liminf / [ f e (u e ) \ d x I(u e ) | 2 dxdt 

JA k Jn 

^ / liminf / f t {u e ) \ d x I{u e ) | 2 dxdt 

JA k e_K) Jn 

^ fc | Afc | . 

So | Ak |< and the set 

£ P;liminf f e {u e ) \ d x I(u e ) | 2 dx = oo| 
has measure zero. We deduce that for all t £ P 

liminf [ I d x I(u e ) I 2 dx < oo 

Jn 

and so for all t £ P 

u e (t ,.) —*• u(t ,.) in H^ +1 (tt) — weakly. 

In particular, we can pass to the limit in the flux J e = f e (u e )d x I(u e ) and write 
lim J e = J = f(u)d x I(u) in L 1 (II) and for almost t £ (0,T). 

e—>-0 

Finally, since u £ u x (t,x ) = 0 for x £ dfi and almost every t £ (0, T). 
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Case 1 ^ a < 2. We recall that (<9 t « e ) is bounded in L 2 (0 So Aubin’s lemma 
implies that ( u e ) converges in L 2 (0,T; C 0,/3 (f2)) for all (3 < 1. We can thus find a subsequence, 
also denoted (it e ), and a set P C (0, T) such that | (0 ,T) \ P | = 0 and for all t £ P , it e (f,.) 
converges strongly in C^{fl). 

We note that for all t £ P, u is strictly positive. Indeed if there exists (to,Xo) S P x such 
that u(to, xo) = 0 then for any (3 < 1 there exists a constant cp such that for all x £ fl 

u(t 0 ,x) < cp | x - x 0 | /3 • 


Thus 


G{u(x , to))dx > 


rdx. 


J J (cp | £ - X 0 \P) n ~ 2 ' 

Given n > 3, we can choose (3 < 1 such that j3(n —2) > 1 . We deduce 

/ G(u(x, to))dx = oo 


which contradicts HMD- 

The rest of the proof is the same as in the first case. 


A Proof of Proposition 14.11 

Our aim is to prove that if 


I us(t,x 2 ) - us(t,x\) |< K | x 2 - Xi | 7 


(45) 


for all t £ (0,T),Xi and X 2 £ fl with constant K independent of <5 and T, then there exists a 
constant M independent of 6 and T such that 

| u s {t 2 ,x 0 ) - usit i,x 0 ) |< M | t 2 - h (46) 


for all ti and t 2 £ (0, T), x £ fl. This proof is an adaptation of the proof done by Bernis-Friedman 
in case 7=^0 Lemma 2.1] for a general 7 . 

We suppose that for all M > 0 one can find xo £ fl and t 2 ,ti £ (0, T) such that 

I u s (t 2 ,x 0 ) - us{ti,x 0 ) |> M | t 2 - ti . (47) 

We suppose that us(t 2 ,x 0 ) > us(ti,Xo) and that t 2 > t±; thus 

us(t 2 ,x 0 ) - us{ti,x 0 ) > M(t 2 - tiY, 0 < fi <t 2 <T, (48) 


where p, = 3 ^+ 3 - We have 





(49) 


where h$ = u^d x I(us), which is valid for any reasonable testfunction. Consider a testfunction Lp 
of the form 


V?(t,x) = Z{x)6 p (t) 
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where £ and Q p are defined as follows. 


£0*0 = £o 


X — Xq 


( M/AK)y (t 2 - h) i 


where M is from (jl51) and K is from (H5l) . and £o0*0 = £o(~ x), £o £ C'o°(fi), £o0*0 = 1 if 
0 < x < £o0*0 = 0 if x > 1 and £' 0 (x) < 0 if x ^ 0. Thus 


£ 0*0 = 


0 if | x — xq |> (M/4K)i (t 2 — ti)T 

k l if | a; - a; 0 K 3 (M/4/v)“ (t 2 -h)^. 


We take 


i if | t-t 2 |<p 


9 (t) = J 0' p (s)ds where 6^(f) = < if | t - h \< p 


0 elsewhere, 


and p < |(t 2 — t\). So, we get 


us£,{x)0'At) = ~ hsi'{x)6 p {t). 


The left-hand side satisfies 

JJ ug(t, x)£(x)8' p (t) -> 4 J £(x)(ug(t 2 , x ) - ug(ti,x))dx as p ->• 0 

To estimate the last expression, we shall only consider values of 2 ; such that 

| x - x 0 |< (M/4:K)~ (t 2 - * 1)2 . 

For such values, 

ug(t 2 ,x) - ug(h,x) = [ug(t 2 ,x) - ug{t 2 ,x 0 )] + [ug(t 2 ,x 0 ) - U,5(tl,Xo)] + [ug(ti,x 0 ) - ug(t 1 ,*)] 

> -2 K \x-x 0 y +M (t 2 - txY 
M, 

> 

Hence, by assuming that the set {£ = 1} is included in Q and by a change of variables in x, 

my 


. 1 . M M a /x 

£(x)(ug(t 2 ,x) - Ug(ti,x))dx ^ / £o(x)dx — (t 2 - - T (t 2 


On the other hand, we have 

hg£'(x)9 p (t) 


€ 


X — Xq 

-T- 


But £'(*) = (. M/AK)y (t 2 - h)~ ) £' 1-r-,, 

in L 2 [Q) we have 


(eo P Y 

so since hg is uniformly bounded 


hgi'{x)9 p {t) 


< 


C 


M A 


(4 K)a 


T {t 2 - ti)y 


2 M 2 ~i 


mi 


— (h - h ) 2 r (t 2 -ti-2p) 2 . 
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Thus by letting p —> 0 we conclude that 


M 1+ 7 (f 2 - < CM~^(t 2 - 

where C is a new constant independent of S, T and M, thus 

M < (t 2 - . 

2~y 

Since fi = 2^+35 we find that M ^ C 3 * 2 ^, and the lemma follows. 
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